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ABSTRACT 


A  new  cross-flow  model  for  a  three-dimensional  momentum  integral 
turbulent  boundary  layer  calculation  method  has  been  developed.  This  cross- 
flow  model  utilizes  two  parameters  and  is  capable  of  representing  the  rever¬ 
sal  (s-shaped)  cross-flow  profiles  that  do  occur  in  ship  boundary  layers. 

The  new  method  uses  only  an  extra  algebraic  equation  for  determining  the 
extra  cross-flow  parameter.  The  basic  momentum  integral  equations  used  by 
the  method  are  the  streamwise  and  cross-flow  momentum  equations  and  the  en¬ 
trainment  equation.  The  method  has  been  programmed  for  calculating  ship 
boundary  layers.  Certain  higher  order  geometric  effects  have  been  included. 
Test  calculations  for  the  experimental  boundary  layer  data  of  the  SSPA-Model 
720  have  been  made.  The  results  are  very  promising.  Good  agreement  with 
experimental  data  is  obtained  over  that  portion  of  the  hull  where  viscous- 
inviscid  interactions  do  not  occur.  Near  the  stern  the  predictions  are 
qualitatively  correct.  The  s-shaped  cross-flow  profiles  are  predicted  near 
the  stern  but  they  are  somewhat  inaccurate.  The  addition  of  viscous-inviscid 
interaction  capabilities  to  the  calculation  should  greatly  improve  the  pre¬ 
dictions  of  the  boundary  layer  near  the  stern. 
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Section  1 


INTRODUCTION 


There  have  been  several  momentum  integral  methods  formulated  for 
calculating  ship  boundary  layers  (cf.  ref.  [l-6]).  All  these  method  have  in 
common  that  the  momentum  integral  equation  along  the  potential  flow  stream¬ 
lines,  the  entrainment  equation  and  the  momentum  integrcl  equation  along  the 
equipotential  lines  on  the  hull  (the  cross-flow  equation)  are  solved  numeri¬ 
cally  for  the  three  main  boundary  layer  parameters.  These  parameters  are  the 
streamline  direction  momentum  thickness  0^,  the  shape  parameter  H  and  the 
angle  3  between  the  surface  shear  stress  vector  and  the  potential  flow  stream¬ 
lines.  The  minor  differences  (although  they  are  important)  between  such 
momentum  integral  methods  consist  mainly  of  the  empirical  formulas  that  are 
used  to  relate  the  boundary  layer  parameters  0^,  H  and  0  to  the  many  other 
parameters  that  are  required  to  describe  completely  the  boundary  layer.  For 
example,  the  surface  skin  friction  coefficient  Cf  is  assumed  to  be  related  to 
the  momentum  thickness  and  shape  factor  H  by  an  empirical  formula.  The 
skin  friction  formula  used  by  various  methods  do  not  differ  in  a  major  way. 
There  are  also  important  differences  among  the  methods  in  the  way  the  equa¬ 
tions  are  solved  numerically.  However,  the  differences  in  the  numerical 
treatment  are  deemed  to  be  technical  and  not  fundamental  ones. 

The  major  fundamental  difference  among  all  the  momentum  integral 
ship  boundary  layer  methods  is  in  the  way  the  cross-flow  is  handled.  The 
standard  way  to  deal  with  the  cross-flow  is  to  assume  a  cross-flow  velocity 
profile  shape  in  terms  of  certain  directly  relevant  boundary  layer  parameters 
such  as  0.  Then  the  main  distinguishing  feature  of  a  method  is  the  nature 
of  the  assumed  cross-flow  velocity  profile. 

This  report  describes  the  results  of  a  research  program  that  had 
several  objectives.  The  first,  most  important  objective  was  to  develop  a  new 
cross-flow  velocity  profile  that  allowed  the  retention  of  the  simplest  form  of 
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the  momentum  integral  formulation  (i.e.,  only  three  differential  equations 
need  be  solved)  but  one  that  encompasses  an  adequate  description  of  s-shaped 
profiles.  Another  objective  was  to  develop  a  numerical  method  to  solve  the 
boundary  layer  equations  which  retained  the  basic  streamline  formulation  but 
did  not  require  a  stream! ine-equi potential  line  grid  system.  Retaining  the 
streamline  formulation  of  the  boundary  layer  equations  is  deemed  important  in 
cases  when  free  surface  effects  must  be  taken  into  account.  In  such  cases, 
the  intersection  of  the  free  surface  and  the  ship  hull  is  a  streamline  which 
forms  one  boundary  of  the  computational  domain.  However,  it  has  been  found  [2] 
that  the  streaml ine-equi potential  line  grid  is  not  a  very  good  system  for  the 
numerical  treatment  of  the  boundary  layer  equations. 

The  third  objective  to  include  higher  order  geometric  terms  in  the 
boundary  layer  equations  is  a  required  prelude  to  the  development  of  a  vis¬ 
cous-  inviscid  interaction  method  for  more  accurate  calculation  of  the  boundary 
layer  close  to  the  stern  of  a  ship. 

Section  2  of  this  report  contains  a  brief  review  of  variors  cross- 
flow  velocity  profiles  that  have  been  developed  in  the  past.  In  Section  3  the 
new  cross-flow  profile  developed  in  this  research  and  its  implementation  is 
discussed.  In  Section  4  the  boundary  layer  equations  including  the  higher 
order  geometric  terms  are  given.  Section  5  discusses  the  new  numerical  method 
for  solving  these  equations.  The  results  of  the  application  of  this  method  to 
the  SSPA  Model  720  [6,7]  test  case  are  described  in  Section  6. 
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Section  2 


» 


CROSS-FLOW  PROFILES 


The  first  moderately  successful  method  for  modeling  the  three-dimen¬ 
sional  turbulent  boundary  layer  was  Mager's  [8].  In  this  momentum  integral 
method  the  cross-flow  boundary  layer  velocity  distribution  v(z),  where  z  is 

*  the  coordinate  normal  to  the  boundary  layer,  was  assumed  to  have  the  form 

jj  =  $  (1  -  |)2tanB  (1) 

*  where  U  is  the  potential  flow  velocity  at  the  edge  of  the  boundary  layer  (here¬ 
after  called  the  slip  velocity),  u(z)  is  the  boundary  layer  velocity  distribu¬ 
tion  in  the  inviscid  streamline  direction  and  6  is  an  overall  boundary  layer 
thickness.  The  velocity  profile  (1)  has  become  known  as  the  Mager  profile  and 

*  it  is  still  the  basis  of  most  of  the  simpler  three-dimensional  turbulent  bound¬ 
ary  layer  calculation  methods.  The  basis  of  this  profile  is  only  three  bound¬ 
ary  conditions.  Profile  (1)  satisfies  the  condition  that  the  cross-flow  veloc¬ 
ity  is  zero  at  the  wall  and  the  edge  z  =  6  and  that  the  limting  streamline 

f  makes  the  angle  8  with  respect  to  the  inviscid  streamline,  i.e., 

limit  77  =  tan&  .  (2) 

z  0  u 

*  The  only  unknown  parameter  of  the  Mager  profile  (1)  is  the  angle  8. 
Hence,  the  single  cross-fl  iw  momentum  equation  is  sufficient  to  determine  8. 

It  has  been  recognized  by  most  researchers  of  three-dimensional 

*  boundary  layers  that  the  profile  (1)  is  unsuitable  for  describing  cross-flow 
velocity  distributions  that  often  occur  in  regions  just  downstream  of  stream¬ 
line  inflection  points.  For  example,  it  has  been  found  often  [4,7]  that  the 
cross-flow  profile  is  s-shaped  in  certain  regions  on  ship  hulls.  The  profile 

*  (1)  is  of  one  sign  across  the  boundary  layer  from  z  =  0  to  z  =  6.  Hence, 
it  cannot  describe  s-shaped  profiles.  In  order  to  encompass  the  variety  of 


» 
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profile  shapes  that  do  occur,  it  seems  evident  that  the  cross-flow  profile 
requires  for  description  more  parameters  than  just  $.  With  the  addition  of 
each  extra  parameter  to  the  cross-flow  profile  an  extra  equation  is  required 
to  determine  this  parameter.  The  addition  of  such  equations  often  complicates 
the  momentum  integral  method  to  the  point  where  the  differential  boundary  layer 
eauations  are  as  easily  solved.  Furthermore,  each  additional  equation  requires 
additional  experimental  data  for  closure. 

There  have  been  many  multi  parameter  cross-flow  profiles  developed. 

The  most  extensive  set  of  these  profiles  was  developed  by  Eichelbrenner  [9]. 

A  similar  set  was  developed  by  Shanebrook  and  Hatch  [10]  .  The  three-dimen¬ 
sional  boundary  layer  methods  based  on  these  profiles  have  not  proved  to  be 
any  more  accurate  than  the  simple  Mager  method  [ll]  .  The  reason  for  this  may 
be  that  these  profile  families  rely  on  parameters  which  are  related  to  the 
slope  of  the  velocity  profile  at  the  outer  edge  of  the  boundary  layer.  It  is 
difficult  to  obtain  reliable  test  data  on  which  to  base  a  differential  eaua- 
tion  that  governs  the  behavior  of  such  parameters. 

Okuno  [5]  developed  a  cross-flow  profile  that  simply  multiplies  the 
Mager  profile  (1)  by  a  linear  function  of  z/6.  This  linear  function  has  one 
free  parameter  which  is  determined  by  the  cross-flow  moment  of  momentum  inte¬ 
gral  equation.  New  unknown  boundary  layer  quantities,  such  as  the  moment  of 
the  shear  stress,  are  introduced  by  this  equation.  Empirical  data  is  lacking 
for  such  terms.  Hence,  the  range  of  usefulness  of  such  an  eouation  is  limited. 
Okuno  does  show  that  his  method  gives  better  predictions  than  Mager' s  profile 
of  the  cross-flow  near  the  stern  of  a  ship  model.  Although  this  method  is  an 
improvement  over  Mager's  method,  the  gain  seems  slight  in  comparison  to  the 
increase  in  complexity  introduced  by  the  cross-flow  moment  of  momentum  integral 
equation. 


It  seems  that  the  construction  of  two  parameter  cross-flow  profiles 
such  as  Eichelbrenner  [9]  and  Okuno's  [5]  do  not  capture  the  essential  physical 
mechanism  that  governs  the  cross -flow  in  the  boundary  layer.  In  the  next  sec¬ 
tion  of  this  report  the  proper  physical  model  of  the  cross-flow  is  described. 


2-2 


I 


t 


A  two  parameter  cross-flow  profile  and  a  simple  algebraic  governing  equation 
for  the  extra  parameter  is  then  synthesized  based  on  this  model. 
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A  NEW  CROSS-FLOW  PROFILE 


The  driving  force  of  the  boundary  layer  cross-flow  is  the  geodesic 
curvature  of  the  inviscid  streamline  on  the  surface  of  the  hull.  This  dimen¬ 
sionless  curvature  term,  denoted  here  by  Kj,  acts  as  a  negative  pressure 
gradient  in  a  direction  to  the  concave  side  of  the  streamline.  Since  the 
velocity  in  the  boundary  layer  is  diminished  from  the  free  stream  velocity, 
the  fluid  particles  there  have  insufficient  momentum  to  withstand  this  pres¬ 
sure  gradient.  Thus,  they  tend  to  flow  down  the  pressure  gradient,  that  is 
towards  the  concave  side  of  the  streamline.  The  fluid  particles  nearest  the 
ship  surface  have  the  least  momentum.  Hence,  a  change  in  sign  of  Kg  causes 
these  particles  to  reverse  rapidly  their  direction  of  motion.  The  fluid  par¬ 
ticles  away  from  the  immediate  vicinity  of  the  ship  surface  are  less  con¬ 
strained.  The  prolonged  application  of  a  pressure  gradient  Kg  to  these  par¬ 
ticles  endows  them  with  sufficient  momentum  that  a  reversal  of  the  sign  of 
l<2  will  result  in  a  belated  reversal  of  the  direction  of  motion  of  these  par¬ 
ticles.  The  disparity  in  the  response  to  changes  of  the  sign  of  the  cross-flow 
pressure  gradient  <2  causes  the  development  of  s-shaped  velocity  profiles  down¬ 
stream  of  streamline  inflection  points  (K2  =  0).  It  is  important  to  note  that 
the  inner  part  of  the  boundary  layer  velocity  responds  quickly  to  the  changes 
in  sign  of  the  pressure  gradient  Kg.  In  fact,  a  careful  observation  of  the 
values  of  B  measured  in  the  experiments  of  Reference  7  and  the  calculated 
values  of  !<2  shows  them  to  be  fairly  closely  correlated.  The  angle  B  is  a 
direct  measure  of  the  direction  of  the  velocity  of  the  inner  boundary  layer 
particles. 


The  cross-flow  momentum  equation  governs  the  momentum  flux  defect 

2 

U  02i  of  the  cross-flow.  The  major  contribution  to  this  momentum  flux  comes 
from  the  outer  part  of  the  boundary  layer.  This  suggests  that  the  cross-flow 
velocity  profile  consists  of  an  extensive  outer  region  which  is  characterized 
by  the  parameter  622*  In  addition  to  this  outer  region  of  the  cross-flow 
there  is  a  small  inner  region  which  is  characterized  by  the  parameter  6. 


Furthermore,  the  cross-flow  velocity  profile  must  satisfy  the  same  boundary 
conditions  that  were  satisfied  by  Mager's  profile  (1).  The  following  profile 
is  suggested: 

S  =  (3) 

where 

*1  u  (1_§) 
o 

and 

r  2 

*2  *  J  U  ^  ^  d^' 

o 

The  cross-flow  velocity  profile  (3)  satisfies  identically  the  definition  of  the 
cross-flow  momentum  thickness 


o 

The  term  proportional  to  in  equation  (3)  gives  the  main  part  of 


the  cross-flow  profile.  The  term  proportional  to  tang  in  equation  (3)  corrects 
the  inner  part  of  the  profile  to  acquire  the  correct  direction  there  (tang). 
Equation  (3)  encompasses  certain  s  shapes  under  the  condition  that  tang  and 
02 j  have  the  same  sign. 

The  cross-flow  momentum  integral  equation  governs  the  distribution  of 
the  quantity  e21  along  the  streamlines  on  the  hull.  It  has  already  been  ex¬ 
plained  that  the  velocity  of  the  fluid  particles  in  the  wall  region  of  the 
boundary  layer  responds  directly  to  the  local  value  of  the  cross-flow  pres¬ 
sure  gradient  l^.  Thus,  it  is  postulated  that  tang  satisfies  the  simple  equa¬ 
tion 


tang  =  cK2 


(1-c) 


e21 

617 


(5) 
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where  c  is  a  coefficient  to  be  determined  by  comparision  of  the  results  using 
this  equation  with  experimental  results.  The  value  of  c=  0  reduces  the  profile 
(3)  to  Mager’s  profile  (1)  and  reduces  the  equation  (5)  to  the  same  relation¬ 
ship  between  tang  and  that  is  obtained  using  Mager's  profile. 

Equation  (5)  properly  accords  with  the  various  ways  in  which  the 
cross-flow  may  occur.  For  example,  a  two-dimensional  boundary  layer  may  be 
suddenly  perturbed  by  a  lateral  motion  of  the  boundary.  Then  02}  is  initially 
zero  and  the  lateral  velocity  of  the  boundary  can  be  related  to  an  effective 
value  of  K2-  In  the  initial  instance  of  such  a  perturbation  the  value  of  tang 
is,  in  fact,  directly  proportional  to  <2.  On  the  other  hand,  the  boundary 
layer  may  develop  first  along  curved  stramlines  and  then  flow  onto  straight 
streamlines  for  which  K2  =  0.  The  boundary  layer  cross-flow  will  then  be 
proportional  to  0^  alone  as  modeled  by  equation  (5).  In  fact,  if  the  stream- 
wise  pressure  gradient  is  adverse  alone  straight  streamlines  then  the  cross- 
flow  thickness  e21  will  grow  in  magnitude  even  with  K2  =  0.  This  condition  is 
encompassed  by  equation  (5).  Note  that  the  opposite  signs  of  021  and  tang 
are  a  consequence  of  the  definition  (4).  An  interesting  phenomenon  is  that 
the  cross-flow  thickness  can  grow  even  after  K2  changes  sign  and  the  inner 
part  of  the  cross-flow  is  driven  in  the  opposite  direction  to  the  outer  part. 
This  may  haopen  in  an  adverse  streamwise  pressure  gradient.  In  fact  this 
occurs  at  the  stern  of  ships  and  is  the  main  cause  of  the  strongly  s-shaped 
cross-flow  velocity  profiles  there. 


9 


Section  4 


THE  BOUNDARY  LAYER  MOMENTUM  INTEGRAL  EQUATIONS 


It  is  not  necessary  to  give  the  detailed  derivation  of  the  boundary 
layer  equations  because  this  derivation  is  tedious  though  elementary.  Instead, 
only  the  final  equations  are  given  and  all  the  terms  are  explained  in  detail. 
The  classical  three-dimensional  momentum  integral  equations  which  do  not  con¬ 
tain  the  higher  order  geometric  terms  can  be  found  in  many  sources  (for  exam¬ 
ple,  references  [1-5]).  The  modified  momentum  integral  equations  that  include 
geometric  correction  factors  will  be  given  in  terms  of  the  thin  boundary  layer 
thicknesses  (such  as  momentum  thickness  and  displacement  thickness)  instead 
of  the  thick  boundary  layer  parameters  such  as  momentum  area,  displacement 
area,  etc.  The  reason  for  this  is  most  of  the  auxiliary  experimental  data  that 
is  needed  to  close  the  system  of  equations  is  given  in  terms  of  the  classical 
thicknesses.  The  conversion  between  the  thicknesses  and  areas  is  entirely 
trivial  as  long  as  the  surface  curvature  is  not  zero.  However,  when  the 
curvature  is  zero,  the  momentum  area,  etc.,  are  undefined  (when  these  are 
originally  defined  in  terms  of  the  local  radius  of  curvature,  as  for  a  body 
of  revolution).  Regions  on  the  surface  for  which  the  total  curvature  is  zero 
are  extensive  on  most  ships;  namely  flat  sides  and  bottoms.  Thus,  it  is  con¬ 
venient  for  ship  applications  to  deal  directly  with  thicknesses. 


Let  s  be  the  arc-length  coordinate  along  an  inviscid  surface  stream¬ 
line.  Let  p  be  the  arc-length  along  a  surface  equipotential  line.  The  momen¬ 
tum  integral  equation  along  the  streamlines  is 
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The  boundary  layer  entrainment  equation  is 
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The  cross-flow  momentum  integral  equation  is 
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TT-  ■  ’•“cre  ‘  2621  ( 1J  § '  “Ki)  -  en(1  +  h)mK2  •  (8> 

The  factors  and  Kj  are  the  geodesic  curvatures  of  the  equipotential  lines 
and  streamlines  respectively.  The  entrainment  parameter  q  is  defined  by 


q  = 


(9) 


where  G  is  Head's  shape  factor.  Head's  shape  factor  G  is  related  to  H  empir¬ 
ically.  The  empirical  data  will  be  given  later.  The  boundary  layer  thicknesses 


are  defined  as 

follows: 
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The  shape  factor  H  is  defined  by  the  equation 


H  =  fij/ejj  .  (11) 

The  terms  and  are  the  surface  skin  friction  coefficients  in  the  stream¬ 
line  and  equipotential  line  directions  respectively.  Equation  (2\  which 
defines  tang,  Implies  that 


Cf2  —  tang  . 


(12) 


I  4 


4-2 


I 


I 


t 


c 


The  function  F(G)  is  the  entrainment  rate.  It  is  given  empirically.  The  terms 
a,  p,  y,  e  and  w  are  the  higher  order  geometric  correction  factors. 


Let  i  be  the  unit  vector  in  the  direction  opposite  to  the  ship  veloc¬ 
ity.  Let  e^  be  the  unit  vector  tangent  to  a  surface  streamline  and  let  Tp  be 
the  unit  vector  tangent  to  a  surface  equi potential  line.  Then  by  definition 
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Kl  =  Ts  *  Tp* 
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K2  =  Vi#- 
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The  geometric  parameters  a,  p,  y,  e  and  w  are  defined  by 
h 

“  '  h  +  t5 
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and  the  factor  h  is  the  normal  radius  of  curvature  of  a  surface  eaui potential 
line.  This  normal  radius  of  curvature  is  defined  by 
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The  terms  t,  u  and  a  in  equations  (13a-e)  are  defined  by 
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The  geometric  terms  a,  p,  y,  e  and  w  that  are  included  in  equations 
( 6 ) - (8)  are  approximations  of  the  most  important  terms  that  occur  in  the  integrals 
of  the  boundary  layer  when  the  variations  of  the  coordinates  in  the  direction 
normal  to  the  surface  are  accounted  for. 

The  following  empirical  data  is  used  to  close  the  system  of  equations 
(6)-(8).  This  data  is  based  on  Head's  ^entrainment  method  for  turbulent  bound¬ 
ary  layers.  The  streamwise  skin  friction  coefficient  is  given  by  the  formula 


Cfl  =  exp(CH  +  B)  (16a) 

C  =  0.019521  +  D(-0. 386768  +  0(0.028345  -  0.000701D))  (16b) 

B  *  0.191511  -  D(0. 83489  +  D(-0. 062588  +  0.001953D))  (16c) 

D  =  ln(RLUen)  (16d) 

The  rate  of  entraimnent  function  F(G)  is  given  by  the  formula 
F(G)  =  0 . 0306 ( G  -  3. Or0,653  (17) 


where 


G  =  1.535(H  -  0.7)“2*715  +  3.3  .  (18) 

The  cross-flow  velocity  profile  has  been  given  by  equation  (3).  The 
coefficient  c  in  equation  (5)  is  given  tentatively  the  value 

c  =  0.125  (19) 

This  value  for  c  seems  to  give  promising  results  in  the  test  calculation  for 
the  SSPA-Model  720.  However,  an  extensive  program  of  calculation  and  comparison 
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with  the  many  sets  of  ship  model  boundary  layer  measurements  that  are  available 
should  be  undertaken  before  a  definitive  value  of  this  coefficient  is  estab¬ 
lished. 


The  cross-flow  boundary  layer  thickness  6^  and  0^  can  be  obtained 
in  terms  of  parameters  e^,  H,g  and  ^  making  use  of  the  cross-flow  veloc¬ 
ity  profile  (3),  the  definitions  (10c,  d)  and  the  following  streamwise  velocity 
profile  u/U  (see  Kool  ^): 

{j  -  [1  -  E(1  -  fAa  -  e‘Qz+)  (20a) 

where 

E  =  1  -  22e‘H/2  v/cTT  (20b) 

(20c) 
(20d) 
(20e) 

n  =  (20f) 

where  L  is  the  scale  length  of  the  ship,  RL  =  (U^D/v  is  the  Reynolds  number, 

Uw  is  the  forward  speed  of  the  ship  and  v  is  the  kinematic  viscosity  of  the 
water. 


Q  =  *sCfl/(l  -  E) 
ZR, 


Z  6 
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Equations  (3)-(20)  establish  a  sufficient  amount  of  information  to 
calculate  the  underwater  boundary  layer  on  a  ship  hull  if  the  ship  hull  geometry 
and  the  Inviscic)  flow  around  the  hull  are  specified.  These  items  are  discussed 
in  the  next  section  on  the  numerical  method  used  to  solve  the  boundary  layer 
equations  (6)-(8).  Certain  terms  involving  the  square  of  the  cross-flow  veloc¬ 
ity  have  been  dropped  from  these  equations  because  they  are  negligibly  small. 

See  references  [l-6]  for  the  complete  boundary  layer  momentum  integral  equations 
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Section  5 


NUMERICAL  METHOD  OF  CALCULATION 


The  numerical  problem  of  calculating  the  boundary  layer  consists  of 
computing  the  inviscid  flow  field  about  the  ship  hull  and  then  numerically 
integrating  the  three  first-order  partial  differential  equations  (6 )- (8)  along 
the  streamlines.  This  report  is  concerned  only  with  the  boundary  layer  problem 
so  it  is  assumed  that  the  inviscid  flow  field  for  zero  Froude  number  flow  (the 
flow  about  the  double  model)  is  known.  It  is  desirable  to  have  a  versatile 
method  for  calculating  the  boundary  layer  that  keeps  the  calculation  of  the 
inviscid  and  viscous  flow  fields  as  independent  as  possible.  Hence  it  is 
assumed  that  the  values  of  the  inviscid  velocity  are  given  at  a  set  of  NP 
points  on  each  of  M  cross-sections  of  the  hull.  The  number  NP  may  vary  from 
one  cross-section  to  another.  The  hull  surface  is  described  by  a  modified  form 
of  the  mapping  method  of  von  Kerczek  and  Tuck^.  Let  x,  y,  z  be  the  coor¬ 
dinates  in  the  longitudinal,  lateral  and  vertical  directions  with  respect  to 
the  centroid  of  the  ship's  design  load  waterplane.  The  positive  x  axis  is  in 
the  direction  of  the  stern  of  the  ship.  The  part  of  the  ship  hull  below  the 
water  is  described  by  the  following  surface  equation: 


R*(x,e)  =  xi  +  y(x,e)T  +  z(x,e)lT 


(21) 


where  (i,j,k)  is  the  set  of  unit  vectors  in  the  x,  y,  z  directions  respectively. 
The  functions  y(x,e)  and  z(x,e)  are  derived  from  the  following  mappings  of  a 
cross-section  of  the  hull  onto  the  lower  half  of  the  unit  circle  in  the  6-plane. 
The  first  mapping  is  the  bilinear  fraction  power  map 


w  -  b  _  K  -  b 
w  +  b  ~  U  +  b 


) 


1/a 


(22) 


where  b  *  b(x)  is  the  starboard  load  waterline  profile,  a=2(l  -  tt/T).  T  is  the 
angle  the  load  waterline  tanqent  of  the  cross-section  makes  with  respect  to  the 
horizontal  y-axis,  c  *  y  +  iz,  i  =  yfT  and  w  is  an  auxiliary  complex  plane.  The 
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purpose  of  this  mapping  is  to  straighten  the  angle  r  to  be  90°  in  the  w-plane. 
Then  the  mapped  cross-section  in  the  w-plane  is  mapped  onto  the  unit  circle  in 
the  f-plane  by  the  method  of  von  Kerczek  and  Tuck^4^  as  follows: 

(23) 

n=l 

where  Re  f  =  0  and  Im  f  =  8.  The  coefficients  an>  n=l,  ....  N,  the  half-beam 
b  and  the  exponent  a  are  all  computed  for  a  number,  say  K,  of  stations  along 
the  hull.  These  values  are  then  interpolated  along  the  hull  using  cubic  Hermite 
polynomial  splines  in  which  the  first  derivative  with  respect  to  x  is  continu¬ 
ous.  This  description  of  the  hull  surface  alleviates  the  two  main  inaccuracies 

fl4*| 

associated  with  the  method  of  von  Kerczek  and  TuckL  J .  The  fractional  mapping 
(22)  allows  any  waterline  slope  instead  of  only  the  vertical  slope  allowed  by 
the  method  of  von  Kerczek  and  Tuck^.  The  cubic  Hermite  spline  lengthwise 
interpolation  of  an's,  a  and  b  considerably  reduces  the  hull  waviness  that 
results  from  the  high  degree  of  polynimial  interpolation  of  the  von  Kerczek 
and  Tuck*-1  method. 

The  availability  of  the  surface  equation  representation  (21)  of  the 
hull  provides  a  natural  intrinsic  surface  coordinate  system.  This  system  is 
made  up  of  longitudinal  lines  of  constant  6  and  the  cross-sections  x  *  constant. 
This  surface  coordinate  system  is  very  nearly  orthogonal  (see  reference  2).  It 
is  visibly  non-orthogonal  only  very  close  to  the  bow  and  stern  of  ordinary 
ships  such  as  the  SSPA  model  720.  On  a  parallel  middle  body  this  surface  coor¬ 
dinate  system  is  orthogonal. 

The  boundary  layer  equations  (6)-(8)  are  solved  on  a  rectangular  grid 
in  the  x-8  plane  that  covers  the  underwater  portion  of  the  hull  surface.  How¬ 
ever,  these  equations  are  integrated  locally  along  the  streamlines.  The  method 
of  this  calculation  can  be  illustrated  by  reference  to  Figure  1  which  shows  the 
rectangular  grid  in  the  x-°  plane 
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Figure  1.  Schematic  of  grid  in  the  x-6  plane  showing  streamline  segments 
along  which  boundary  layer  equations  are  integrated. 

The  boundary  layer  calculation  procedure  begins  at  a  station  xQ, 
designated  as  k  *  kQ.  There  the  primary  boundary  layer  parameters  e^,  q  and 
e21  are  given  (say  by  experiments)  on  a  set  of  points  6j.  The  given  values  of 
©22 »  9  and  e21  are  su^icient  t0  calculate  the  values  of  all  the  other  required 
boundary  layer  parameters  using  the  various  formulas  of  Section  4. 

The  inviscid  flow  velocity  is  known  at  the  cross-section  k  on  some 

o 

other  set  of  points.  These  values  of  the  inviscid  velocity  are  first  inter¬ 
polated  with  respect  to  the  coordinate  e,  using  cubic  Hermite  splines,  onto  the 
points  0..  These  values  of  the  inviscid  velocity  are  then  used  to  calculate 

J 

the  segments  of  the  streamlines  that  pass  through  the  points  0j  at  station  kQ 
and  intersect  the  cross-section  at  kQ  +  1  at  the  points  designated  by  the  6j* 
as  shown  in  Figure  1.  The  inviscid  velocity  at  station  k0  +  1  is  then  inter¬ 
polated  to  the  streamline  intersections.  Equations  (6)-(8)  can  now  be  inte¬ 
grated  along  the  streamline  segments  from  the  points  at  station  kQ  to  the 
respective  streamline  intersections  at  the  station  kQ  +  1.  The  integration 
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method  that  has  been  used  for  this  purpose  is  one  that  evaluates  the  cross¬ 
derivatives  3S2/3p  and  3612/3p  at  station  kQ.  Then  equations  (6)-(8)  are  inte¬ 
grated  forward  one  step  as  ordinary  differential  equations,  with  respect  to 
s,  using  a  perdictor  corrector  method  based  on  the  Euler  and  Trapezoidal  inte- 
gragion  rules  (see,  for  example,  reference  15).  After  the  values  of  6^,  q 
and  have  been  calculated  at  the  streamline  intersection  points  at  station 
kQ  +  1,  these  values  are  interpolated,  again  using  cubic  Hermite  splines,  onto 
the  points  ft..  The  procedure  can  then  be  carried  out  exactly  in  the  same  way 

J 

to  integrate  the  equations  (6)-(8)  from  station  k0  +  1  to  kQ  +  2,  and  so  on 
to  the  end  of  the  ship. 

This  method  has  proved  to  be  very  easy  to  implement  and  has  the  fol¬ 
lowing  convenient  features.  The  number  and  distribution  of  integration  points 
ft.  at  each  station  can  be  changed  easily  at  every  step.  This  can  accommodate 

J 

regions  of  rapid  changes  of  the  boundary  layer  parameters  and  allows  marching 
around  local  regions  of  separation.  Either  or  both  of  the  side  boundaries  (in 
the  example  shown  in  Figure  1,  the  keel  and  load  waterline)  can  be  boundary 
streamlines  and  need  not  follow  a  line  of  constant  ft.  For  example,  if  the  com¬ 
plete  ship  problem  is  treated  in  which  the  free  surface  is  not  assumed  to  be 
flat,  then  the  load  waterline  streamline  is  replaced  by  the  free  surface  stream¬ 
line.  The  boundary  layer  calculation  method  does  not  have  to  be  carried  out  on 
the  same  set  of  points  as  the  ones  at  which  the  inviscid  velocity  is  sDecified. 
The  only  requirement  is  that  the  inviscid  velocity  is  specified  at  some  set  of 
points  on  each  cross-section  k.  This  can  be  accomplished  easily  by  a  suitable 
interpolation  preprocessing  of  whatever  inviscid  velocity  data  is  provided  (for 
example,  by  a  Hess-Smith  type  of  calculation). 

The  above  description  of  the  calculation  procedure  focused  directly 
on  the  main  elements  of  the  numerical  method  for  integrating  the  equations  (6)- 
(8)  along  the  surface  of  the  hull.  No  mention  was  made  of  the  calculation  of 
all  the  geometric  data  such  as  Kj,  K^,  h,  dp,  ds,  directional  unit  vectors  and 
inviscid  velocity  gradients  in  order  to  keep  the  discussion  as  simple  as  possi¬ 
ble.  The  required  geometric  and  pressure  gradient  data  is  fairly  easy  to  cal¬ 
culate  using  the  surface  equation  (21).  For  example,  the  unit  vector  e*n  normal 
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to  the  surface  is  calculated  directly  by  analytically  differentiating  the  sur¬ 
face  equation  (21).  The  unit  vector  Ts  tangent  to  a  streamline  is  easily  cal¬ 
culated  using  the  inviscid  velocity.  Then  the  unit  vector  e*p  along  the  surface 
equi potential  lines  is  calculated  by 


ep  =  enxes- 

The  arclength  increments  ds  and  dp  are  easily  determined  using  the  surface 
equation  (21)  and  the  vectors  e*s  and“ep.  The  the  derivatives  3U/9s,  de^/ds 
and  dejj/dp  are  calculated  by  difference  approximations.  These  calculations  are 
very  similar  to  the  ones  described  in  reference  2.  The  reader  is  referred  to 
this  reference  for  the  details. 

The  important  aspect  of  this  boundary  layer  calculation  method  that 
makes  possible  a  fairly  easy  way  of  calculating  all  the  required  surface  geom¬ 
etry  is  the  representation  of  the  ship  hull  surface  by  the  conformal  mapping 
functions.  This  method  of  representing  the  hull  surface  has  been  discussed  in 
detail  in  reference  14.  However,  the  modification  using  the  fraction  power 
mapping  (22)  is  new  so  that  it  is  worthwhile  to  discuss  briefly  why  this  map¬ 
ping  is  used  and  the  improvements  that  are  gained  in  the  surface  representation. 
This  can  be  illustrated  best  by  referring  to  Figure  2.  This  figure  shows  the 
given  input  offsets  for  station  x  =  +0.9L,  where  L  is  the  half-length  between 
perpendiculars,  of  the  SSPA-Model  720  and  the  representation  of  this  station 
by  the  conformal  mapping  method  with  and  without  the  preliminary  mapping  (22). 
The  mapping  (23)  used  N  =  7  coefficients  in  both  cases.  The  vertical  slope  at 
the  load  waterline  (z  =  0)  of  the  section  representation  without  the  mapping 
(22)  cannot  be  eliminated  by  choosing  larger  values  of  N.  It  is  easy  to  see 
that  by  using  the  preliminary  mapping  (22)  the  representation  of  the  section 
is  dramatically  improved. 

This  numerical  method  has  been  implemented  in  the  computer  program. 
The  detailed  description  of  this  implementation  is  left  for  an  appropriate 
Program  Manual.  The  next  section  describes  the  results  of  a  test  calculation 
for  the  SSPA-Model  720. 
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Figure  2.  Comparison  of  the  mapping  representation  of  section  x  =  0.9L  of 
the  SSPA-Model  720  with  input  offsets.  The  input  offsets  are 

the  circles.  The  dashed  curve  ( - )  is  the  representation  by 

the  mapping  formula  (23)  alone.  The  solid  curve  ( - )  is  the 

representation  by  the  combined  mapping  formulas  (22)  and  (23). 
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COMPUTATIONAL  RESULTS 


The  boundary  layer  calculation  method  described  in  Section  3  through 
5  was  programmed  and  tested  on  the  SSPA-Model  720.  This  is  a  double  model 
of  a  moderate  block  coefficient  ( Cg  =  0.675)  cargo  hull.  Extensive  boundary 
layer  measurements  were  made  on  this  model  by  Larsson^.  This  model  was 
used  as  one  of  the  test  cases  to  which  many  different  boundary  layer  calcu¬ 
lation  methods  have  been  applied  in  the  SSPA-ITTC  Workshop  on  Ship  Boundary 
Layers^.  Reference  to  the  Workshop  report^  can  be  made  for  all  the  details 
of  the  model  and  the  experimental  data,  as  well  as  the  comparison  of  the  per¬ 
formance  of  various  boundary  layer  calculation  methods. 

The  computed  boundary  layer  by  the  present  method  will  be  compared 
with  the  experimental  data  for  the  SSPA-Model  720.  No  direct  comparisons  of 
the  present  method  with  any  of  the  other  methods  presented  at  the  Workshop 
will  be  made  because  this  would  require  too  much  replotting  of  the  published 
data.  Reference  to  the  report ^  can  be  made  to  compare  the  present  method 
with  the  others  of  the  Workshop. 

The  present  boundary  layer  calculation  method  requires,  like  all 
other  methods,  the  specification  of  certain  boundary  layer  parameters  along  an 
initial  cross-section  of  the  hull.  The  boundary  layer  can  then  be  computed 
for  stations  downstream  of  the  starting  cross-section.  The  present  calcula¬ 
tion  method  requires  the  specification  of  four  boundary  layer  parameters  at 
the  starting  station.  The  most  convenient  parameters  to  be  specified  are  the 
streamwise  momentum  thickness  0^,  the  shape  factor  H,  the  cross-flow  angle  B 
and  the  cross-flow  momentum  thickness  02j.  However,  starting  data  at  station 
x  =  -0.5L,  where  L  is  the  half-length  between  perpendiculars  and  the  origin 
of  x  is  at  midship,  given  in  reference  7  is  the  distribution  of  6^,  H  and  B 
only.  Hence  the  starting  distribution,  of  the  values  of  Sjj  were  estimated  on 
the  basis  of  the  relationship 
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(see  equation  5).  Equation  (24)  is  the  result  of  assuming  that  the  Mager 
cross-flow  velocity  profile  (1)  is  valid.  This  approximation  for  the  initial 
distribution  of  e 21  may  he  responsible  for  part  of  the  discrepancy  between 
the  computed  and  the  measured  boundary  layer  characteristics  discussed  below. 

The  boundary  layer  downstream  of  station  x  =  -0.5L  was  computed  for 
a  dense  set  of  points  on  the  hull  surface  (see  Figure  1).  It  is  simplest  to 
discuss  the  results  by  comparing  them  with  measured  data  along  certain  curves 
on  the  surface  of  the  hull.  The  first  set  of  comparisons  between  computed  and 
experimental  results  is  the  distrubition  of  values  of  ©11  and  B  along  two 
streamlines,  designated  streamlines  B  and  C  (see  reference  7).  Streamline  B 
runs  down  the  hull  close  to  the  keel  except  near  the  stern  where  it  turns 
upward  across  the  bilge  towards  the  location  where  the  propeller  shaft  would 

emerge  from  the  hull.  Streamline  C  runs  down  the  hull  on  the  center  of  the 

bilge  and  emerges  at  the  stern  front  the  hollow  beneath  the  stern  section  flare. 
Figure  3  gives  a  comparison  between  the  distribution  of  computed  and  experi¬ 
mental  values' of  the  streamwir;  momentum  thickness  6^  along  streamlines  B 
and  C.  This  figure  also  shows  the  experimental  values  of  611  along  the  keel, 
designated  as  streamline  A,  for  comparison  with  the  values  of  e along  stream¬ 
line  B.  Direct  calculation  of  the  boundary  layer  along  the  keel  was  not  made. 
It  can  be  seen  in  Figure  3  that  the  values  of  6^  are  well  predicted  by  the 

present  method  on  streamline  C,  but  on  streamline  B  the  values  of  6^  are 

considerably  overpredicted.  Calculations  of  the  boundary  layer  were  also  made 
using  a  small  cross-flow  approximation  to  equations  (6)-(8).  In  this  approxi¬ 
mation  the  terms  involving  6^  and  6^  in  equations  (6)  and  (7) ,  respectively, 
are  dropped.  The  cross-flow  has  no  effect  on  the  streamwise  flow  in  this 
approximation.  The  resulting  distribution  of  values  of  6^  along  streamline 
B  is  not  substantially  different  from  the  values  shown  in  Figure  3.  Hence 
the  overprediction  of  the  values  of  0^  on  streamline  B  cannot  be  attributed 
to  cross-flow  effects.  It  seems  likely  that  the  values  of  the  pressure  gradi¬ 
ents  along  streamline  B  calculated  by  the  present  method  may  be  larger  than 
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the  actual  values  in  the  experiment.  The  results  shown  in  Figure  3  are  very 
encouraging  though. 


Figure  3.  Comparison  of  the  experimental  and  computed  values  of  the  stream- 
wise  momentum  thickness  along  streamlines  B  and  C  on  the  SSPA- 
Model  720.  The  curves  are  the  computed  values  and  the  points  are 
the  experimental  values  of  en. 


Figure  4  gives  a  comparison  of  the  distribution  of  computed  and  ex¬ 
perimental  values  of  the  cross-flow  angle  0  along  streamlines  B  and  C.  Here 
the  similarity  between  the  computed  and  experimental  results  is  heartening. 
The  discrepancy  between  the  computed  and  measured  values  of  0  near  the  start¬ 
ing  station  x  =  -0.5L  on  streamline  B  may  be  due  to  the  possibly  inaccurate 
initial  values  of  02}  estimated  by  using  formula  (24).  It  should  be  noted 
that  most  of  the  values  of  0  are  fairly  small.  The  measured  values  of  0  were 
obtained  by  measuring  the  inner  part  of  the  boundary  layer  velocity  profiles 
using  a  Preston  tube^.  The  calibration  and  even  the  applicability  of  this 
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technique  for  measuring  the  flow  direction  close  to  the  wall  is  uncertain. 
Hence,  a  precise  similitude  between  the  computed  and  measured  values  of 
should  not  be  expected. 


Figure  4.  Comparison  of  the  experimental  and  computed  values  of  the  cross- 
flow  angle  8  along  the  streamlines  B  and  C  on  the  SSPA-Model  720. 
The  curves  are  the  computed  values  and  the  points  are  the  experi¬ 
mental  values. 

Figures  5,  6  and  7  show  the  distribution  of  values  of  6^  and  B 
on  the  cross-sections  at  x  =  0.5L,  0.71  and  0.9L,  respectively.  The  computed 
and  measured  values  of  9^  and  8  are  in  satisfactory  aoreement  at  station 
x  =  0.5L.  At  station  x  =  0.7L  the  comparison  between  the  computed  and  meas¬ 
ured  values  is  less  satisfactory  but  not  too  bad.  At  station  x  =  0.9L  the 
computed  values  of  6^  and  8  are  only  qualitatively  similar  to  the  respective 
measured  values.  The  main  part  of  the  discrepancies  between  computed  and 
measured  values  of  6^  and  6  shown  in  Figures  6  and  7  is  most  likely  due  to 
the  occurrence  of  viscous-inviscid  interactions  that  have  not  been  accounted 
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•  for  in  the  present  computational  method.  A  second  factor  causing  this  dis¬ 

crepancy  may  be  the  very  poor  resolution  of  the  given  pressure  distribution 
near  the  keel  at  station  x  =  0.9L.  The  interpolation  near  the  keel  of  the 
potential  flow  data  given  by  the  Workshop  is  very  erratic.  This  causes  the 
f.  values  of  0  to  be  negative  and  very  large  in  magnitude  near  the  keel  at  sta¬ 

tion  x  =  0.9L. 

0,,  x  103 
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Figure  5.  Comparison  of  the  distribution  along  the  cross-section  of  the  com¬ 
puted  and  experimental  values  of  the  streamwise  momentum  thickness 
011  and  the  cross-flow  angle  0.  In  this  figure  the  cross-section 
is  at  x  =  0.5L. 
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Figure  6.  Comparison  of  the  distribution  along  the  cross-section  of  the  com¬ 
puted  and  experimental  values  of  the  streamwise  momentum  thickness 
9ll  and  the  cross-flow  angle  0.  The  cross-section  for  this  data  is 


» 


Figure  7.  Comparison  of  the  distribution  along  the  cross-section  of  the  com¬ 
puted  and  experimental  values  of  the  streamwise  momentum  thickness 
6ll  and  the  cross-flow  angle  0.  The  cross-section  for  this  data 
is  x  *  0.9l- 
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It  is  interesting  to  note  that  only  one  method  tested  in  the  SSPA- 
ITTC  Workshop  seemed  to  give  substantially  better  results  at  station  x  =  0.9L 
than  the  present  method.  This  was  the  differential  method  of  Soejima, 

Yamazaki  and  fakatake  (see  reference  7).  This  method  made  use  of  higher  order 
boundary  layer  theory  in  which  the  viscous-inviscid  interactions  were  accounted 
for.  The  results  shown  in  Figures  6  and  7  do  indicate  that  with  the  addition 
of  viscous-inviscid  interaction  capabilities  the  present  method  may  be  capable 
of  predicting  well  the  boundary  layer  characteristics  near  the  stern  of  a 
ship. 

Viscous-inviscid  interaction  effects  would  tend  to  smooth  the  vari¬ 
ation  of  the  pressure  distribution  across  the  cross-sections.  Thus  the  dis¬ 
tribution  of  values  of  the  boundary  layer  parameters  would  also  be  less  vari¬ 
able  across  the  cross-sections.  This  would  be  in  accord  with  the  experimental 
data  shown  in  Figures  6  and  7. 

Figures  8  and  9  show  the  profiles  for  the  streamwise  and  cross-flow 
velocity  at  station  x  =  0.9L  on  streamline  B  and  C.  The  computed  profiles  are 
not  very  accurate,  but  are  at  least  in  qualitative  agreement  with  the  experi¬ 
mental  profiles.  In  particular,  the  cross-flow  Drofile  on  streamline  B  is  of 
the  reversing  (s  shape)  type.  The  lack  of  quantitative  aareement  between  the 
computed  and  experimental  velocity  profiles  should  not  be  attributed  to  the 
formulas  for  the  profile  shapes.  This  lack  of  agreement  is  simply  due  to 
the  fact  that  the  boundary  layer  parameters  9^,  H  and  6  were  not  predicted 
accurately  enough  at  station  x  =  0.9L.  It  is  believed  that  inclusion  of  the 
effects  of  viscous-inviscid  interactions,  not  any  basic  change  in  the  boundary 
layer  cross-flow  model,  will  bring  the  computed  values  of  the  boundary  layer 
parameters  into  agreement  with  experimental  values. 
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Figure  8.  Comparison  of  the  computed  and  experimental  velocity  profile 
on  streamline  B  at  station  x  =  0.9L.  u/U  is  the  velocity  in 
the  streamline  direction  and  v/U  is  the  velocity  in  the  cross 
flow  direction.  The  curves  are  the  computed  profiles. 
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Figure  9.  Comparison  of  the  computed  and  experimental  velocity  profile  on 
streamline  C  at  station  x  =  0.9L.  u/U  is  the  velocity  in  the 
streamline  direction  and  v/U  is  the  velocity  in  the  cross-flow 
direction.  The  curves  are  the  computed  profiles. 
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CONCLUDING  REMARKS 
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This  report  describes  the  development  of  a  new  three-dimensional 
momentum  integral  method  of  calculating  ship  boundary  layers.  It  has  been 
shown  that  the  method  does  give  reasonably  good  predictions  of  the  boundary 
layer  on  the  part  of  the  ship  where  thin  boundary  layer  theory  is  valid. 

The  results  of  the  test  calculation  near  the  stern,  while  not  unreasonable, 
do  require  improvement.  Such  improvements  will  result  by  including  viscous- 
inviscid  interaction  effects.  The  basic  boundary  layer  model  seems  adequate. 

The  main  accomplishment  of  this  effort  has  been  a  fairly  simple 
generalization  of  the  Maper  cross-flow  model  to  include  s-shaped  profiles. 
What  is  needed  now  is  to  validate  this  model  by  further  calculation  of  well- 
documented  test  cases.  Such  a  program  of  calculations  win  help  to  establish 
the  value  of  the  constant  (or  possibly  the  function)  c  in  equation  (5)  that 
relates  the  cross-flow  angle  8  to  the  cross-flow  pressure  gradient  <2  and  the 
momentum  thickness  62 j. 
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